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^ ! Abstract 

in 



We obtain the hole propagator of the Sutherland model with SU(2) internal 
symmetry for coupling parameter (3=1, which is the simplest nontrivial case. 
One created hole with spin down breaks into two quasiholes with spin down 
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and one quasihole with spin up. While these elementary excitations are ener- 
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getically free, the form factor reflects their anyonic character. The expression 

C*~> ' for arbitrary integer is conjectured. 
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There has been a remarkable development in the study of one dimensional models of 
particles with the inverse-square interactions |I]. Stimulated by the pioneering work by 
Simons, Lee, and Al'tshuler ||, the dynamical correlation functions of the Sutherland model 
for spinless particles were obtained explicitly with the help of various techniques, such as 
the supermatrix method ||, 2D Yang-Mills theory ||, and the theory of the Jack symmetric 
polynomial [HH. 

The amazing success in the dynamics is, however, limited to the spinless case and much 
less is known for that of the multicomponent system. The multicomponent Sutherland model 
(MSM) [^,|| realizes the multicomponent Tomonaga-Luttinger liquid in the simplest man- 
ner. The deeper understanding of the MSM would be helpful to understand other strongly 
correlated electron systems, especially t-J model [H and the singlet factional quantum Hall 



effect flIO | . Thus the study on the dynamics of the MSM is highly desirable. 

In this Letter, we present the results on the one-particle Green function of the Sutherland 
model for particles with SU(2) spin (<Tj = ±1). This model is the simplest example of 
the MSM. This is the first calculation of the dynamical correlation function of integrable 
continuum models with internal symmetry. 

We consider the following Hamiltonian: |J7| 

- f <9 2 27T 2 P(P + P tJ ) 

t.dxy I? fZstffrixi-xJ/L] { ] 

for iV-particle system. Here P^ is the operator that exchanges the spin of particles % and 
j fLl| |. The dimensionless coupling parameter (3 takes non- negative integer values. L is the 



linear dimension of the system. The statistics of particles are chosen as boson (fermion) for 
odd (even) (3, so that we can set Py = (— 1) M^- with the coordinate exchange operator 

Ma H- 

For the model ([j]), we calculate the hole propagator: 

G? (x, t) = (g, JV|$ (x, t) ^ T (0, 0) |g, N)/(g, N\g, N) (2) 

for (3 = 1, which is the simplest nontrivial case. Here ip-^(x,t) represents e lWi ?/>f(x)e _lWt . 



|g, N) is the ground state for iV-particle system. The main result of this paper is the 
following explicit expression in the thermodynamic limit with p = N/L: 

G**{x,t) = ai f 1 dn f 1 dv 2 f 1 dv ^-^ /3 ^-^r 2/3 \^-^r 2/3 

U J-i J-i [ ei K) ei (t; 2 ) ei (*; T )] 1/3 

x exp [-it {ei (vi) + e t (v 2 ) + e 1 (/u t ) - Ci} + i (vi + v 2 + u-f) vxpo/2] , (3) 

with 01 = poCiF [2/3] /(72r 2 [1/3]), Ci = (Snp /2) 2 , and e, (v) = 3 (np /2) 2 (1 - v 2 ) . 

In the following paragraphs, we present the derivation of @. In the last part of this 
paper, we conjecture the expression for arbitrary non-negative integer (3 on the basis of the 
result for (3 = 1. We consider only the case where N is even and M = N/2 is odd. 

The ground state is singlet and its wavefunction $^ ({X, = exp[i27r£j/L], crj}) is given 

by0] 

$f ({X t , a,}) = Uf <3 {(*< - X j f +S{ ' Ti '^ (i)^"^)/ 2 } 

,^^-iHM-^,,^^ (4) 

s 

with Pg [{<7i}} = Hies $ [fh 1] Tljgs $ [ a jy ~~ !]■ The sum with respect to 5 runs over the subset 

of {1, 2, 3 • • • , N} with M elements. The norm of (f§ is given by N C M D(M 1} M 2 ; [3)L N , with 

N C M = N\/[M\(N - M)\] and 

d (m u M 2; # = n / p- n / 4^- n i* - ^r 2 n k - ^r 2 n n i* - ^-r . 

(5) 
Here C is the contour counterclockwise on the unit circle. The value of D (Mi, M 2 ; (3) was 



derived in ref. |T2 



as 



with 



^"^n^ff, w 



A(J 1 ,J a ;/?) = (J a l)n r[(/3 + l) J + J 1 + l] fo^^^" 1 - ( 7 ) 



When J 2 = 1, the product is taken as unity. 



We consider the action of if) on the ground state. The wavefunction of the state 
V>t(0,0)|g,iV) is given by 

N 



'£(-l)W>*>[({X i ,<r i },X k =l,* k = l). (8) 

fc=i 

Since particles are identical, we pick out the term k = N in Eq. (g): 

JV-l 

$f ({X„ a,,} ; X* = 1, a N = 1) = (~i) M II ^ (*« - l)*^"' 11 ^~\ (9) 



t=i 



where $^ _1 represents the ground state wavefunction with momentum — n(M — 1)/L. We 
calculate the propagator for Eq. (|]) and then multiply it by N. The factor N comes from 
the number of ways to choose which particle is removed and added. 

For taking the trace over spin variables in Eq. (|2]), we replace the Hamiltonian ([[]) by 

* ti dxl I?h sin 2 [tt ( Xi - x 3 ) /L] ■ [ U) 

Spin variable for each particle is then a good quantum number since Hamiltonian ( |10D does 
not contain spin variables. Hence we can fix spin for each particle and multiply the result by 
n-\Cm-\ as an equivalent procedure to trace over spin variables. It is from this advantage 
that we consider ( |i0|) instead of ([I]). In the following, we set <7j = 1 for 1 < i < M — 1 and 
Oi = -1 for M < i < 2M - 1. 

After tracing over {&i}, we arrive at the following expression: 

JV-l / „ iv \ Af-1 M 



G'ix, t) = co n [f ^-J *r n & - ^ +1 n ^ - v" 

x e -i{(«,-^)«-P*} J| (^ _ !)/»+! -Q (w . _ 1} /3 ^ (n) 



with Co =jv_i C m-iPo/[nCmD (Ml, M2,/3)]. The overlines represent complex conjugate. 
Here $^ f ~ 1 is given by 

Af-1 M Af-1 M 

^r 1 = n (* - ^ +1 n k - ^f +l n n c* - ^ 

i<jf fc</ i=l jf=l 

Af-1 Af 

X II ^(^H^/ 2 JJ w -[/3(AT-l)+M-l]/2 ) (12) 

i=l i=l 



where Zi (1 < i < M — 1) = X t and Wj (1 < j < M) = Xj + m-i represent the complex 
coordinates of particles with spin up and down, respectively. 

Instead of calculating G f3=l (x, t) directly, we first calculate the following correlation func- 
tion for later convenience: 

7V-1 / „ jv- \ Af-1 M 



&{x, t) = co n (/ ^-j ^r 1 n & - i) 2 n & - 1) 

x e-^-^y-^ n c* - 1) 2 n (^ - *) ^ r ~ 1 - ( i3 ) 

t=i i=i 

From (0) to (0), we just replace fli Ilj 0* - l)"" 1 " 1 Oj - I)' 3 by n» Ilj (^ - !) 2 (™j ~ 1) but 
leave the rest of (3 as an arbitrary integer. 

The first step in calculating &(x,t) is to expand Ufii 1 (zi ~ tflljLi (wj ~ 1) $f _1 in 
terms of eigenf unctions of the Hamiltonian Tip. Relevant eigenstates have the following 
form: 

K{{z ilWj })^-\ (14) 

Here K({zi,Wj}) represents a polynomial satisfying the following conditions: 

• symmetric with respect to exchange between Zj's and exchange between w/s 

• polynomial of degree not more than 2 with respect to Zi and of degree not more than 
1 with respect to Wj. 

Now we assign two types of indices to eigenstates described by (|14]). The polynomial part 
K ({zi,Wj}) is written by the linear combination of rh^ ({zi,Wj}) = m^ ({<^}) m.„j. ({ w j})- 
Here m^ ({zi}) and ra^i ({wj}) are monomial symmetric functions of {zi} i= ^ and {wj},_ v 
respectively. These are given by 



A/-1 T M i 



"fy 



Pi P j 



respectively. The sum with respect to P runs over distinct permutations. Here \x is defined by 
H = {yu^/i 1 }, where /i T = {/*L A*I» A*s ' ' - A*L-i} and ^ = {/*x» A*a» /*i — , A*iwr} are partitions 



arranged in decreasing order {ft[ > ft > ft 3 - ■ ■ , ft > ft > ft • ■ •). The weight J2fLi 1 ft 
(EjLiHj) is denoted by |/x T | (\ft\). 

Next we define ft and ft by ft = ft + {M — 2, M — 3, • • • , 0} and ft = /i^ + 
{M — 1, M — 2, • • • , 0}, respectively. Furthermore we introduce /i + as the rearranged se- 
ries of ft © ft in decreasing order. Here we define the order of ft' so that fi + > v + if the 
first nonvanishing difference fxf — vf is positive. Next we define the order of [i = I ft; ft >. 
If /i + > v + , then fj, > v. If fi + = v + and \ft\ > \v^\, then fi > v. 

We define the polynomial K^ ({zi,Wj}) as 

• Kp ({z h Wj}) = m^ ({z h Wj}) + Eu«ri v^rhv {{z h Wj}) 

• K M ({zi, Wj}) ^g' 1 is an eigenfunction of Tip. 

• K^ ({zi, Wj}) $^ _1 's are orthogonal with one another. 

Eigenf unctions of Tip are obtained for distinguishable particles by the same method as used 
for the single component bosonic case |jl3| ,|14|j . 



Although the index /i seems natural to describe the polynomial K^, here we introduce 
another index A = {Ai,A2,Aj}, which is related directly with rapidities or velocities of 
elementary excitations. Ai (A2) is the smaller (larger ) non-negative integer among the set 
{M, M - 1, • ■ • , 0} - z> T and A T is among {M, M - 1, ■ ■ ■ , 0} - iA. In Fig. 1, we illustrate 
the case of /i = I ft; ft\ = {2211; 1111} for example. We see that there are three vacant 
rapidities A = {M — 5, M — 1, M — 4}. Now we identify Ai,A 2 and A^ as rapidities of 
quasiholes with spin down and up, respectively. 

In order to calculate G 13 , we have to know the expansion coefficient C\ which appears in 



M-l M 

\2 



n {^ - 1) 2 n K- - 1) = E E E c x k x (i6) 

i=\ j=l 0<Ai<A 2 1<A 2 <M0<A T <M 



and the norm N\ for each wavefunction: 



N- 1 r AY 9 

2=1 t 



In terms of C\ and N\, the hole propagator is written as 



&(x,t) = c J2 C'^exp [-i (E*- 1 - E» ) t + iP; 



\x 



(18) 



Ai,A 2 ,Af 



Now we present the explicit forms of E^ 1 , P\, C\, and N\. Momentum of the state A is 
given by 



P A = 7T (5M -l)/L- (2tt/L) (Ai + A 2 + A T ) + tt/3 (1 - 2M) /L. 



In terms of A = {Ai, A2, A-]-}, we write the expression P A — P^ as 



(19) 



E N-1 _ E N 



(2tt/L) 2 [- (M - 1/2) Mf3 2 - (M 2 - 10M + l) /4 — /3 (2M 2 + 7M - l) /2 + (2M - 1) /? 2 /4 
- (1 + 2/?) (A 2 + A 2 + A 2 ) + (M + 2/3M + p - 1) (Ai + A 2 + A T ) - 2(3 (% - A 3 )] . (20) 

Here Ai, A 2 , and A 3 are the smallest, the second smallest and the largest of {Ai,A 2 ,A|}, 
respectively. 

General expressions of C\ and N\ for the Jack polynomial have been derived in math- 



ematical literature | I5| . In the case of K\, however, these expressions have not yet been 
derived. Hence we calculate K\ for system with finite number of particles [M < 11) up 
to degree |//| = 20. Based on these data, we obtain the general form of C\ and N\. The 
coefficient C\ is given by 



c A = r 2 [(i + /3)77]/(r[(2 + /3)77]r 

' V[ x 2 -x 1+ ri\ r[A T -A!+(l-/3)t?] r[A T -A 2 +r?] 
r[A 2 -Ai-/3?7] r[A T -A!+r;J r[A T -A 2 +(l+/3)??J 

for Ai < A 2 < Af 

r[A 2 -Ai+(l-/3)t?] r[A T -Ai+?y] t[\ 2 -\ ] -2(3tJ\ 
r[A 2 -Ai-2,3j7] r[A T -Ai + (l+/3)r ? ] r[A 2 -A T -/37j] 

for Ai < At < A 2 

r[A 2 ^A 1+1? l r[Ai-A T -2/3T?] r[A 2 -A T -3^ 



x < 



r[A 2 -Ai-/3r;] r[Ai-Ai~ 0rj\ T 

for At < Ai < A2, 



A 2 -A T -2/3?7 



(21) 



with rj = 1/(1 + 2/5). The norm is given by 



N x 



D(M - 1, M; j3)T [M] T [M + 1 - /fy] .A A g A , 

r [m - /?r/] r [m + 1 - 3/377] r 3 [(i + /?) ^ ' 



(22) 



with 



3 r 



Afc + 1 — k(3r) 



kii r [A fc + 1 + (1 - fc)/fy 



M-A fc + (l + (A;-2)/3)77 



M - A fe + (1 + (A; - 1)/3)t7 



(23) 



and 



#a = 



r[A 2 -A 1 -/3> ? ]r[A 2 -A 1 +(l+/3)>?]r 2 [A T -A 1 +l-2/3> ? ]r 2 [A T -A 2 +l-/3r;] 

r[A2-Ai]r[A 2 -Ai+r ? ]r[A T -Ai+l-3/3» ? ]r[A T -Ai+l-/3r;]r[A T -A2+l-2/3r;]r[A T -A2+l] 



for Ai < A2 < At 



r[A 2 -A 1 +»?]r[A 2 -Ai-2/3r ? ]r 2 [A T -A 1 + (l+/3)r ? ]r 2 [A 2 -A T -/3r;] 

r[A 2 -Ai+(l-/3)» 7 ]r[A 2 -Ai- / 3» ? ]r|A T -Ai+l|r|A T -Ai+» ? |r|A 2 -A T -2/3r,|r|A 2 -A T | 



for Ai < At < A 2 



r[A 2 -Ai-/3> ? ]r[A 2 -Ai+(l+/3) } 7]r 2 [Ai-A T -/3r ? ]r 2 [A 2 -Aj-2/3> ? ] 
r[A 2 -Ai]r[A 2 -Ai+r?]r[Ai-A T Jr[A 1 -A T -2/3» ? Jr[A 2 -A T -/3» ? Jr[A 2 -A T -3 / 3r,J OT T < * < 2 ' 



(24) 



Now we consider the thermodynamic limit M — ► 00, L — *■ 00 with p — 2M/L. For 
Ai, A2, A-f, |Ai — A2I, |Ai — A-f], I A2 — A|| ~ C(M), we obtain the asymptotic forms of C\, 
N\, E^~ l — Eg , and Pa- In the case of C\ and N\, we can use the Stirling formula: 
r [N + 1] — »■ y/2nN N+1 / 2 e~ N . The sums with respect to Ai, A2, and A-f turn into the multiple 
integral: 



2 \ 3 A 2 -l MM 1 fl /"I /■! 

with Vi — 1 — (2Aj/M) (2 = 1, 2, f ). From the above results, &(x, t) reduces to 

1 1 11 |2(l+/3)r) I |-2/?n i i — 2/3n 

G^,t) = a,/ 1 d, 1 / 1 d, 2 / 1 d,^ U21 ^- Vl1 '^-^ 



(25) 



[e/3 (vi) e/3 (u 2 ) e^ (v T )] /3?? 
x exp [-it { e/3 ( Vl ) + e p (v 2 ) + ep (v T ) - (,3} + i (2 - 2(3 + Vl + v 2 + v T ) 7nrp /2] , (26) 



with the prefactor 



Q/3 



2 i+6^ ( 2 + pf r 2 [(2 + (3) 7]] r 2 fo] ' 



(27) 



the chemical potential (p = [ir(l + 2j3)p /2} 2 , the quasihole energy ep (t>) = 
(7rp /2) 2 (1 + 2/3) (1 - v 2 ), and the exponent 7 = (5 + (3)r] - (1 + (3). By setting (3 = 1 
in Eq. (pop, we obtain the expression (|3]). 

Now we discuss the physical implication of (p6|). We can see that the three quasiholes 
are energetically free while the form factor has nontrivial structure. From the study of 
thermodynamics [|l^], on the other hand, we have learned that quasiholes are free but obey 
nontrivial (exclusion) statistics. Here we attribute the nontrivial form factor to the exclusion 
statistics. Namely we read the form factor in (p6|) as 

\vi — U2I 9u \v-\ — Vi\ 9n \v-\ — V2I 9n 



[^w^wr^wi 



ii— srt 



(28) 



where the matrix 



9h 



( \ 

9tt 9n 



1 + 2/3 



/ 



(3+1 -13 
-(3 (3 + 1 



(29) 



y 9 it 9u J 

represents the statistical interaction among quasiholes of the Sutherland model with SU(2) 
internal symmetry. Actually we meet the matrix (|29|) in the thermodynamics of the present 
model. 

Lastly we consider the hole propagator G"(x,t) with arbitrary integer coupling (3. In 
this case, we know that one created hole with spin down breaks into (3 + 1 quasiholes with 
spin down and (3 quasiholes with spin up. Energies and momenta of intermediate states 
are obtained easily. In the thermodynamic limit, we obtain E^ 1 — E^ — > J2i=i e p ( v i) + 
EjLi € ( Uj ) - ((3 and Pf- 1 -► (vrpo/2) (E^l v t + £? =1 %) , where Vi (1 < i < (3 + 1) and 
Uj (1 < j ; < /3) are normalized velocities of quasiholes with spin down and up, respectively. 
Although we do not have the expression for the form factor, we make a conjecture relying 
on the interpretation (|2"8]). The resultant expression for G^(x,t) is given by 



/3+1 .1 P 

&{x,t)~i[ / d^n 

x exp 



duj 



n 



/3+1 



i<j \ V i V 3 



.\ 2 9unP 



n^K-^i^n^n 



2<?n n^ +1 n? 



j — l \Vi Uj 



.\ 2 9U 



n&foMf^nUfoM] 1 ' 8 " 






[ 1=1 



i=i 



i=l 



J=l 



(30) 



up to a constant factor. 

The expression ( |30"D would be plausible in view of the known form factor in the single 
component model (with integer coupling (3), which is given by 



IT 9 I _ \ 29 

lU<j \Vi v j\ /oi\ 

nfK^)] 1 " 9 ' [ } 

Here V{ is the velocity and e(v) is the energy of quasihole. The statistical parameter g 
was derived as 1/(3 from the study of the thermodynamics fll7 |. We can see that the form 



factor in ( |3~0"1) is a natural generalization of (0). The conjecture ( ^0|) should be checked 
microscopically, e.g. with the use of the theory of the nonsymmetric Jack polynomial |18 . 
This is left as a future work. 

The author thanks Y. Kuramoto and T. Yamamoto for useful comments. This work was 
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Ministry of Education, Science, Sports and Culture of Japan. 
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FIGURES 
FIG. 1. A diagram that contributes to &. This figure corresponds to the case of 
fi = {2211; 1111}. The open and shaded blocks represent the ground state and an excitated 
state, respectively. There are two quasiholes with spin down at Ai and A2, while one with spin up 
is located at Af. 
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